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Abstract 

In this paper, we prove a local null controllability result for the three-dimensional 
Navier-Stokes equations on a (smooth) bounded domain of R 3 with null Dirichlet bound- 
ary conditions. The control is distributed into an (arbitrarily small) open subset and has 
two vanishing components. J.-L. Lions and E. Zuazua proved that the linearized system is 
not necessarily approximately controllable even if the control is distributed on the entire 
domain, hence it cannot be null controllable and the standard linearization method fails. 
We use the return method together with a new algebraic method inspired by the works of 
M. Gromov and previous results by M. Gueye. 
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1 Introduction 

1.1 Notations and statement of the theorem 

Let T > 0, let be a nonempty bounded domain of K 3 of class C°° and let lj be a nonempty 
open subset of il. We define Q C R x I 3 by 

Q := (0,T) x n = {(t,x)\te (0,T) and x G 0}. 

The current point x € M 3 is x — (xi, x 2 , x 3 ). The i-th component of a vector (or a vector field) 
/ is denoted f l . The control is u = (u 1 , it 2 ,u 3 ) <E L 2 (Q) 3 . We require that the support of u is 
included in to, which is our control domain. We impose that two components of u vanish, for 
example the first two: 

u 1 = and u 2 = in Q, (1.1) 

1 Work supported by the ERC advanced grant 266907 (CPDENL) of the 7th Research Framework Programme 
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so that u will be written under the form (0,0, \ u v) with v G L 2 (Q) from now on, where l u : 
fl R is the characteristic function of to: 

l w = 1 in ui, l w = in O \ u. 

Let us define 

V :={yeHl(nf\V • y = 0}. 

The space 1/ is equipped with the 7J 1 -norm. Let us denote by H the closure of V in L 2 (fl) 3 . 
The space ff is equipped with the L 2 -norm. 

We are interested in the following Navier-Stokes control system: 

Ay + (y ■ V)y + Vp = (0, 0, l u v) in Q, 

V • y = in Q, (1.2) 

y = on[0,T]x9Q. 

From [27, Theorem 3.1, p. 282], we have the following existence result: For every y° g H, there 
exist y £ L 2 {{0, T),V) n £°°((0, T), ff) and p € L 2 (Q) satisfying 

y(0,-)=j,°infi (1.3) 



Vt 



such that (1.2 1 holds. 

Our main result is the following theorem, which expresses the small-time local null-controllability 
of ( |L2| : 

Theorem 1. for every T > and for every r > 0, there exists r\ > smc/i i/iat, /or every y° G V" 
verifying \\y Wh 1 !^) 3 ^ ??, there exist v £ L 2 (Q) and a solution (y,p) G £ 2 ((0, T), i/ 2 (f2) 3 n V) l~l 
L oc ((0,T),fl' 1 (f2) 3 ny) x L 2 (Q) o/ (fL2|)-([Ll|) suc/i ffctrf 

y(r,-)=0, (1.4) 

lklk 2 (Q) 3 «S ^ (1-5) 

H2/||L 2 ((o,T),ff 2 (n)3)nL°°((o,T),ffi(^) 3 ) ^ r - (1-6) 
Remark 1. Once a control v G L 2 (Q) is given, the corresponding solution (y,p) of (1.2), ( |1.3| 
and ( |1.6| given by Theorem^is unique (recall that for the Navier-Stokes system, the uniqueness 
of(y,p) means that y is unique and p is unique up to a constant depending on the time). This 
comes from the uniqueness result given in ]2Tl Theorem 3.4, P- 297]: One has 

^((o.rj.ir^n) 3 ) c ^((cr^fr^fi) 3 ) c l 8 ((o,t),l 4 (^) 3 ) 



thanks to a classical Sobolev embedding, and there is at most one solution (y,p) of (1.2 1 and 
( 1.3 I in the space 

L 2 ((0, T), V) n L°°((0, T),H)n L 8 ((0, T), L 4 (fi) 3 ) x L 2 (Q). 

Remark 2. One observes that in Theorem^the initial condition y° is more regular than usual 
(y G H). In fact, using the same arguments as in ]15] and f!6f (see also |6j Remark 1]), 
one can easily extend the previous theorem to small initial data in H f] L 4 (tt) 3 with a solution 
(y,p) G L 2 ((0,T),V) n L°°((0,T),H) x L 2 {Q). In this case, Rem ark [7| is no longer true and 
there might possibly exist many solutions (y,p) verifying (1.2 I and (1.3| once v is given. 
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1.2 Some previous results 

The controllability of the two or three-dimensional Navier-Stokes equations with a dis- 
tributed control has been studied in numerous papers. In general, for Navier-Stokes equations, 
it is relevant to consider the approximate controllability, the null controllability or the control- 
lability to the trajectories, the second one being a particular case of the third one. 

In [21] , a first result of local exact controllability to the trajectories was established under 
technical conditions: had to be homeomorphic to a ball, the control had to be supported 
in a nonempty open subset whose closure is included in f2, and the target trajectory had to 
be a stationary solution of the Navier-Stokes equation. Moreover, there were some technical 
regularity conditions for these stationary solutions. A similar result for the linearized Navier- 
Stokes equations was established but with the same strong conditions. Many of these hypotheses 
were removed in |22| . 

Then, it was proved in jTS] the local exact controllability to the trajectories with regularity 
conditions that were weaker and more suitable for the study of the Navier-Stokes equations. 
In this article, the authors also proved some exact controllability results for linearized Navier- 
Stokes systems, with very weak regularity conditions. The same authors proved in the 
local exact controllability with a control having one vanishing component, provided that ui 
"touches" the boundary of the domain in some sense. Later on it was proved in |13| a local 
null controllability result on the Stokes system with a control having a vanishing component 
without the geometrical condition on u, but the authors were not able to extend it to the 
nonlinear Navier-Stokes system. A recent work ([5]) improved the previous one and proved 
the local null controllability of the Stokes system with an additional source member by means 
of a control having a vanishing component, which enabled the authors to prove the local null 
controllability of the Navier-Stokes system for a control having a vanishing component. In all 
these articles, the main points of the proof were to establish first the controllability of the 
linearized control system around the target trajectory thanks to Carleman estimates on the 
adjoint of the linearized equation, and then to use an inverse mapping theorem or a fixed-point 
theorem to deal with the nonlinear system. 

The natural question is then: Can we remove another component of the control, which 
would provide an optimal result with respect to the number of controls? Reducing the number 
of components of the control is important for applications, and have already been studied many 
times for linear or parabolic systems of second order (that are quite similar to linearized Navier- 
Stokes systems), see for example [I], where a necessary and sufficient condition to control a 
system of coupled parabolic equations with constant coefficients and with less controls than 
equations is given, or (3j [25] for time-dependant coefficients. If the coefficients depend on the 
time and the space, there are no general results, notably if we consider two coupled parabolic 
systems where the coupling region and the control region do not intersect (a partial result 
under GCC is given in [T]). For a recent survey on the controllability of coupled linear parabolic 
equations, see [5]. 
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1.3 The linear test 

To obtain Theorem ^ the first natural idea is to linearize the system around 0, i.e. to 
consider the Stokes control system 

y t - Ay + Vp= (0,0, ul,,) in Q, 

v • y = in <3. (1 7) 

y(0,-)=y° mil, 
y = on[0,T]xdQ. 

It is well known (see for example [TB], or (22| ) that if this linear system were null controllable 
(with, in addition, an arbitrary source term in a suitable space), then applying an inverse 
mapping theorem (for example the one presented in [2]) in some relevant weighted spaces, 



we would obtain that (1.2 1 is locally null controllable around 0. However the linear control 



system (1.7 1 is in general not null controllable, and not even approximately controllable: In 



it is proved that if we choose Q to be a cylinder with a circular generating set and with 



an axis parallel to e 3 then (1.7) is not approximately controllable, even if we control on the 
entire cylinder. However, it is proved in |24| that the approximate controllability property holds 
generically with respect to the generating set of the cylinder. 

Since linearizing around is not relevant, we are going to use the return method, which 



consists in linearizing system (1.2 1 around a particular trajectory (y,p, u) (that we construct 
explicitly), verifying 1/(0, •) = y(T, ■) — 0, proving that the linearized system (with a source 
term / verifying an exponential decrease condition at time t = T) is null controllable, and then 
concluding by a usual inverse mapping argument that our system is locally null controllable. This 
method was introduced in [7] for a stabilization problem and first used in [8 for controllability 
issues for partial differential equations. The return method was already successfully used in 
HU HHt to obtain global controllability results for the Navier-Stokes equations and in [T5] to 
prove the small-time local null controllability for the Navier-Stokes equations on the torus T2 
when the control has one vanishing component. For more explanations about the return method 
and other examples, see |10l Chapter 6]. 

2 Constructing a relevant trajectory 

In this subsection, we construct explicit particular trajectories (y,p,u) going from to 
so that, as it will be shown in section [3] the linearized control system around them is null 
controllable. 

Without loss of generality we may assume that £ cj. Let g £ C°°(IR 3 ), 

g : (t,w,x 3 ) 1 V g(t,w,x 3 ), 

and h £ C°°(IR 3 ), 

h : (i, w, X3) i-> h(t, w, X3). 
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For (xi,x 2 ,x 3 ) £ M 3 , let r := y/x\ + x\. We define y £ C°°(R 4 ;]R 3 ) by 



y(t,x) 



fg(t,r 2 ,x 3 )x 1 
g(t,r 2 ,x 3 )x 2 I ,Vt € M, V.t = (xx,x 2 ,x 3 ) £ 



\ h(t,r 2 ,x 3 ) 

Let 7"! > be small enough so that 

Ci := {(x 1 ,x 2 ,x 3 ) elV^ ri, |x 3 | < n} C w. 

On the functions g and /i, we also require that 

Supp(. g ) c [T/4,T] x (-00, r?] x [-n,n], 
Supp(fc) c [T/4,T] x (-oo, r?] x [-n,n]. 



(2.1) 

(2.2) 

(2.3) 
(2.4) 



In (2.3 



(2.4 1 and in the following, Supp(J) denotes the support of the function /. From (2.1 ), 

(2.5) 



(|2.2|), (|2.3j and ( |2.4| , one obtains 

Supp(y) C [T/4, T] x C x C (0, T] x oj c (0, T] x ft 



Let p G C°°(]R 3 ) be defined by 



P(t,w,x 3 ) 



1 



(g t - {iwg ww + 8g w + g 33 ) 



2wgg w + \g\ 2 + hg 3 ) (t, w' , x 3 )dw' . (2.6) 



In ( |2.6[ ), in order to lighten the notations, for I £ W. N £ W and i = (ii, . . . ,ijv), with 
ife G {t,w, 3} for every fc G {1, ...AT}, and 93 : D(<p) C M 3 -> JR J , (t, 10,2:3) i-> ip(t,w,x 3 ), 
</ ? iii 2 ---ijv denotes the following derivative: 



• • • 



where 



£j = x 3 if j = 3, £j = t if j = < and £j = w if j = w. 
We use similar notations throughout the paper. Let p £ C°°(]R 4 ) be defined by 

p(t, X!,x 2 ,x 3 ) := p(t, r 2 , x 3 ). 



From (2.21, (2.3 1, Q2.4|, \2.Q\ and (2.7), it follows that 



Supp(p) c [T/4,T] x d c [T/4, I*] x w c (0,T] x O. 



(2.7) 



(2.8) 
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From pTE) , ( [276] ) and ( |277| , one obtains 

y\- + {y-V)t +n=^ 

f t -Ay 2 + (y-V)y 2 +p 2 =Q. 

Let u G C°°(IR 4 ) 3 be defined by 

5:= (0,0,y?~Ay 3 + (y.V)y 3 +p 3 ). 
From ( |2.11[ ), one obtains fl7l| . From ( |2~5| . ( [278] ) and ( |2.11| , we have 



(2.9) 
(2.10) 

(2.11) 



Supp(u) c (0,T] x w. 



From (2.91, (2.101, (2.11) and (2.121, we have 



Finally, in order to have 



it suffices to impose 



Vt - &y + (y • v )y + = \ u u. 



div y = 0, 



dh 

dx 3 



~2(g + wg w ). 



(2.12) 



(2.13) 



(2.14) 



(2.15) 



Let v be a positive numerical constant which will be chosen later on. Let a G C°°(ffi), 6 G C°°(R) 
and c G C°°(K) be such that 



We choose 



and 



Supp(a) C [T/4,T] and a(t) = e< r -'> 5 in [T/2,T], 
Supp(fe) C (-co, r\) and b(w) = w, Vs G (— 00,7^/4], 
Supp(c) C (-ri,ri) and c(x 3 ) = xj in [-ri/2, r x /2]. 

g(t,w,x 3 ) = sa(t)b{w)c! (x 3 ) 

h(t, w, x 3 ) = — 2ea(i)(£>(u>) + w6'(u>))c(:E3), 



(2.16) 
(2.17) 
(2.18) 

(2.19) 

(2.20) 



where e > will be chosen later on. From ( 2.19[ ) and (2.201, one obtains (2.151. 

In the next section, we prove that, for every small enough T, for every small enough e > and 
for a well-chosen the linearized control system around the trajectory (y,p,u) is controllable. 
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3 A controllability result on the linearized system 

3.1 Definitions and notations 

The linearized control system around the trajectory (y,p,u) is the linear control system 

' y\ - Ay 1 + (y ■ V)?/ 1 + (y ■ V)f + Pl = f in Q, 



y 2 t - Ay 2 + (y • V)y 2 + {y • V)y 2 + p 2 = f in Q, 

y 3 - Ay 3 + (y • V)?/ + {y • V)y 3 + p 3 = l u v + f in Q, 

V • y = in Q, 

y = on[0,T]x<9O, 



(3.1) 



where the state is y : Q — > M 3 , f : Q — >• M 3 is a source term (it will be specified later in which 
space exactly it shall be) and the control is v : Q — > M. Let uj n be a nonempty open subset of 



C 2 := {(x 1 ,X2,x 3 );r < —, \x 3 \ < — | , 
which will be chosen more precisely in the next section. Let 

Qo := (T/2.T) x w - 

The following figure represents the open subsets of ft we introduced up to now. 
Let £: C°°(Q ) 5 -> C°°(Q ) 4 be defined by 

^ y t 1 -Ai/ 1 + (y-V)i/ 1 + (yV)y 1 +pi ^ 
yf - Ay 2 + {y ■ V)y 2 + (y ■ V)y 2 + p 2 
yf - Ay 3 + (y ■ V)y 3 + (y ■ V)y 3 + p 3 - v 
V-y 



(3.2) 




(3.3) 



for every y = (y 1 ^ 2 , y 3 ) G C°°(Q ) 3 , for every p E C°°(Qo) and for every v e C°°(Qo)- Let us 
denote by 

£ := (ar ,a;i,a;2,a;3) = (t,x 1 ,x 2 ,x 3 ) 
the current point in Q - For a = (ao, a\, a 2 , a 3 ) G N 4 and ip : Q — > K fe , 9"</3, denotes, as usual, 

(in this part, this notation will be more convenient than the one introduced in the previous 
section, which is more compact only for derivatives of small order). Let £(M fe ;M z ) be the set of 
linear maps from R k into R l and Mk,i(R-) be the set of matrices of size k x I with value in the 
ring 1Z. 

As usual, in the inequalities written in this article C denotes a constant (depending in general 
only on u>, 0, T) that may change from one line to another. 
Let us give some other definitions. 
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Figure 1: The open subsets Ci,C2,Uq,u. 



Definition 1. A linear map A4 : C°°(Qo) k —> C°°(Qo) 1 is called a linear partial differential 
operator of order m if, for every a — (ctQ, a%, ct2, CK3) £ N 4 with \a\ := ao + ot\ + 012 + «3 ^5 m, 
there exists A a G C°°(Q ; C(R k ;R 1 )) such that 

(M(p)(0= E ^(05>(0,V^eQo,V^eC 00 (Q ) fc . 

I a 1 m 

j4 linear map Ai : C°°((3o) fe — > C°°(Qo)' is called a linear partial differential operator i/ 
i/iere exists m G N smc/i £/ia£ .M is a linear partial differential operator of order m. 

Let fcbea positive integer and let B := (B 1 ,B 2 ,B' i ) : C°°(Q a ) k -> C°°(Q ) 3 be a linear- 
partial differential operator. Let us consider the linear equation 





B 2 f 




B 3 f 


w 


V / 



(3.4) 
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where the data is f E C°°(Qo) k and the unknown is (y,p, v) E C°°(Qo) 5 - Following [19, Section 
2.3.8, p. 148], we adopt the following definition. 



Definition 2. The linear equation (3.4) is algebraically solvable if there exists a linear partial 
differential operator M : C°°(Qo) k — > C°°(Q ) 5 such that, for every f E C°°{Qo) k , Mf is a 
solution of ( 3.4 1 , i.e. such that 



CoM = (B,0). 



(3.5) 



In the following, every function <p E C°°{Qof with a compact support included in Qo is 
extended by in Q \ Qo and we still denote this extension by ip. 

The next proposition explains how the notion of algebraic solvability can be useful to reduce 
the number of controls as soon as a controllability result is already known for a large number of 
controls. In fact, the question of the null-controllability of (3.1) can be split into two distinct 
problems: One "algebraic" part (solving system) 3.5)) and one "analytic" part (finding controls 
which are in the image of B, the control acting possibly on all the equations and not only the 
third one) . This proposition has a very general scope and could be formulated for more general 
control systems. It is inspired by techniques used in the control of ordinary differential equations 
(see, in particular, |10l Chapter 1, p. 13-15]). 



Proposition 1. Let us consider the linear control system 



yf - Ay 3 



(F-V)y 1 
(F-V)y 2 
(F-V)y 3 



(yV)y 2 + P2 =B 2 u 
(y V)y 3 +p 3 



V-y 

y 

2/(0,0 



B 3 u- 



y° 



f 1 
f 
f 



in Q, 
in Q, 
in Q, 
in Q, 

on [o,t] x an, 

in 



(3.6) 



where the state is y : Q — » IR 3 , the control is u E C ca (Q) k , which is required to have a support 
in Qo, and f := (Z 1 ,/ 2 ,/ 3 ) E C°°(Q) is a source term. Let us assume that: 



A\. The linear control system (3.6 I is null controllable during the interval of time [0,T] in the 
sense that for every y° E V and for every f E C°°(Q) such that 



there exists S > such that f = on [T — S,T] x Q, 



(3.7) 



there exists u E C°°(Q) k with a compact support included in Qo such that the solution 
(y,p) of (3.6) with initial condition y(0, •) = y° satisfies y(T, •) = 0. 



A 2 . 



(3.4 1 is algebraically solvable. 



9 



Then, the linear control system (3.1 1 is null controllable during the interval of time [0, T]: For 
every y° G V and for every f £ C°°(Q) satisfying (3.7), there exists v £ C°°(Q) with a compact 



support included in Qq such that the the solution (y,p) of (3.1 1 satisfying the initial condition 
2/(0, •) = y° satisfies y(T, •) = 0. 



Proof of Proposition [TJ First of all, we use the null-controllability of (3.6 I with controls in 
the image of B (Assumption Ai): let y° £ V and let u* £ C°°(Q) k with support included in Qo 



such that the solution {y*,p*) of the Cauchy problem (3.6) for the initial condition y*(0, •) = y 
satisfies y*(T, •) = 0. Let us remark that B is a local operator, which implies that Bu^ c still 



has a compact support in Q . Now we use the algebraic solvability of (3.4 1 (Assumption ,4 2 ) : 
Let M. be as in Definition [5J For a map h £ C°°(Q) k with a support included in Qo, we denote 
by M. the map from Q into E 5 defined by 

Mh = in Q \ Q ,Mh = M(h lQo ) in Q . 

We shall use this slight abuse of notation until the end of the paper. Note that, for every 
h £ C°°(Q) k with a support included in Qo, Mh £ C°°(Q) 5 and has a support included in Qo 
(because M. is a local operator) . Let us call 

{y,Pi v) := —Aiu*. 

One observes that the support of (y,p,v) is still included in Qo- Li particular y(T, •) = 0. Let 

(y,p) ■= (y* +y,p* +P)- 



Then 



y\ - Ay 1 + {y ■ V)^ 1 + (y ■ V)y x + Pl 



Vt 



Ay 2 
Ay 3 



(y V)y 2 +p 2 
(yV)f+P3 
V -y 

y 

2/(0,-) 

y(T,-) 



f 1 
f 

l«w + f 





y° 
o 



in Q, 
in Q, 
in Q, 
in Q, 
on [0, T] x dfl, 
in Q, 
in f2. 



This shows that the linear control system (3.1 1 is indeed null controllable during the interval of 
time [0,T] and concludes the proof of Proposition [l] ■ 



Remark 3. For the sake of simplicity, we have formulated Proposition^ in a C°° setting. Let 
us assume that the control u coming from Assumption A\ is not of class C°° , but is less regular 
(one observes that the regularities of y* , p* and f do not matter for the proof of Proposition^ 
since only u* is differentiated by the linear partial differential operator M ). For example, assume 
that u* £ Hi where Hi is a functional space (for example a weighted Sobolev space), and assume 
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that M can be extended on Hi, Aiu* being then in another functional space Hi (for example 
another weighted Sobolev space of order less that Hi in order to take into account that M. is 
a linear partial differential operator). Then one easily verifies that Proposition^ remains true 
as soon as every function of Hi (and its derivatives until the order at least the order of Ai) 
vanishes at time t = T, the first assumption Ai being changed as the following: the linear control 
system ( |3.6[ ) is null controllable during the interval of time [0,T], i.e. for every y° 6 V and for 
every f e L 2 (Q) satisfying dO) th ere exists u € Hi with support included in Qq such that the 
solution (y,p) of (3.6 I satisfying the initial condition y(0, •) = y satisfies y(T, •) = 0. Similarly, 
we will need to relax property (3.7 1 by replacing it with a suitable decay rate near t = T. This 



will be detailed in Subsection \3.4\ 

It remains to deal, for a suitable choice of B, with Assumption Ai (we shall do it in Subsec- 
tion 3.2 1 and with Assumption Ai (we shall do it in Subsection 3.3 1. 



3.2 Algebraic solvability of (3.4) 
We choose k = 7 and define B by 



B(f\f 2 J 3 J 4 J 5 ,fj 7 ) := 



ffl + fi + fi 
ff + fi 
V f 7 



fi 



The main result of this subsection is the following proposition. 



(3.8) 



Proposition 2. There exists e* > 0, there exists T* > such that, for every e £ (0,£*) ; there 
exists a nonempty open subset loq of such that Assumption Ai holds for every T < T* : There 
exists a linear partial differential operator A4 : C°°(Qo) 7 — ^ C°°(Qo) 5 such that (|3.5[) holds. 



3.2.1 The adjoint problem 

Let £q : C°°(Qo) 4 — > C ca (Qo) 3 be the linear partial differential operator defined by 




hi - Ay 1 + (y ■ V)y x + (y ■ V)y x 
y\ - Ay 2 + (y ■ V)y 2 + (y ■ V)f 




(3.9) 



for every y = (y\y 2 , y 3 ) € C°°{Q ) 3 , and every p € C°°(Q ). 

Let Bq : C°°(Qo) 6 — > C°°(Qo) 3 be the linear partial differential operator defined by 



Bo(/\/ 2 ,/ 3 ,/ 4 ,/ 5 ,/ 6 ) := 



ffl + fi + fi 
ft + fi + f! 
V o 



(3.10) 
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Note that the third equation of (3.4 1 can be read as 

v = y 3 - Ay 3 + (y • V)y 3 + (y • V)y 3 + p 3 - f . 

Hence, one easily sees that Assumption A2 is equivalent to the existence of a linear partial 
differential operator M : C°°(Q ) 6 ^(Qo) 4 such that 



£ oM =B . 



(3.11) 



As in [19, p. 157], we study (3.111 by looking at the "adjoint equation". For every linear 
partial differential operator M. : C°°(Qo) fe ~ > C°°(Qq) l , M. = Yl\ a \<m A a d a , we associate its 
(formal) adjoint 

M* -.C^iQo) 1 ^C°°(Q ) k 

defined by 



|q.| 



(3.12) 



where is the transpose of the matrix A a (£). (Definition (3.12 1 makes sense since J2\ a \<m A a d a — 

implies that the A a are all equal to 0.) One has M** = M and, if M : C co (Q ) k -> C°°(Qo) 1 
and J\f : C°°((3o) Z C°°((5o) m are two linear partial differential operators, then (J\f o M)* = 
M*oM*. 



Hence, (3.11) is equivalent to 



M* oQ=B* . 



(3.13) 



Direct computations, together with (2.14), show that, for every z — (z 1 , z 2 ) £ C co (Qo) 2 and for 
every tt € C°°(Q ), 



a 



(-zt-Az 1 -(y-V)z 1 +y{ 
-z 2 - Az 2 - (y ■ V)z 2 



viz 2 



v\z l + y\z 2 - 7T 2 



Viz 1 + y\z 2 - 7T 3 



V 



(3.14) 



Bl 



. 7T . 



( z \i z 2i Z \, Z 2 , Z 2 7 Z 2 ). 



(3.15) 



Assumption A2 is now equivalent to the following property: There exists a linear partial dif- 
ferential operator N(= M* ) : C°°(<2o) 4 -> C°°(Qo) 6 such that for every tp = {(p 1 ,^ 2 , ip 3 , tp 4 ) € 
C°°(Q ) 4 , if (^ 1 ,z 2 ,tt) e C°°(Q ) 3 is a solution of 



then (z^ , ^2 1 ^3 > %i 1 %2 1 z 3 ) 



-z t x - Az 1 - (y ■ V)z 1 + y\z l + y\z 2 - 7Ti = ip 
-z 2 - Az 2 - (y ■ V)z 2 + y\z X + V\z 2 ~ ^2 = V 

Viz 1 + Viz 2 - 7T 3 = Lp 3 , 

1 2 4 

- Zf = 99 4 , 



(3.16) 
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Remark 4. The most natural linear partial differential operator B to try first would have been 
B : C°°(Q ) 3 -> C°°(Q ) 3 defined by 

f f i> 

Bf:= I f 2 I ,V/=(/ 1 ,/ 2 ,/ 3 )GC 00 (Qo) 3 . (3.17) 



3 



Unfortunately, Proposition^ does not hold with this B. Indeed, in this case Bq : C° c (Qq) 3 



C°°(Qo) 2 would be now (compare with (3.15),) such that, for every z = (z^z 2 ) G C°°(Qq) 2 and 
for every it 6 C°°(Qo), 

B* Q (z,n) = (z\z 2 ). (3.18) 

Let Fx e C°°(T/2,T) and let F 2 G C°°(T/2,T). We define z = (z\z 2 ) G C^IQq-.R 2 ) and 
tt G C°°(Qo) by 

z\t,x) -Frit), z 2 (t,x) :=F 2 {t), 

and 

n(t,x) := -F[{t) Xl - F^(t)x 2 + F^y 1 + F 2 (t)y 2 . 



Then Cq(z,tt) = 0. However, if (i*\,_F 2 ) ^ (0,0), then Bq(z, tt) ^ 0. Hence, in this case, (3.131 
does not hold whatever the linear partial differential operator M.q is and whatever the trajectory 
(y,p,u) is. 

3.2.2 Number of variables and equations 

Let us give some algebraic results about the number of derivatives of a certain order. 

Definition 3. Consider a scalar PDF with a smooth (enough) variable z depending on 4 vari- 
ables Xo, x±, X2, X3. We call equations of level n all the different equations we obtain by dif- 
ferentiating the PDF with respect to all the possible multi-integers of length n. The number of 
distinct equations of level n is denoted E{n), and the number of distinct equations of a level less 
than or equal to n is denoted F(n). 

Remark 5. Clearly, E(n) is also the distinct number of derivatives of order n for (smooth 
enough) functions having 4 variables, and F(n) is also the distinct number of derivatives of an 
order less than n for (smooth enough) functions having 4 variables. Moreover, if we consider a 
scalar PDF with many variables z 1 , . . . , z k depending on xq,xi,X2, £3 containing derivatives of 
z 1 , . . . z k of order m at most, the maximum number of derivatives of z 1 , . . . z k we may expect in 
the equations of a level less than or equal to n is kF(n + m). 

We want to compute E and F precisely. One has 

£(n)= (n + l)(,+ 2 )( , + 3) ; (3ig) 

F(n)= (n + l) ( n + 2)(n + 3)(n + 4)_ ^ 
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Indeed 



(a , a 1 , a 2 , a 3 ) H> {a + 1, ckq + a x + 2, ao + a i + a 2 + 3, a + a 1 + a 2 + a 3 + 4} 



defines a bijection between the set of (ao, ai, ot-z, a.3) £ N such that «o + a i + ^2 + <^3 ^ n 
and the set of subsets of {1, 2, . . . , n + 4} having 4 elements. Hence, .F(n) being the number of 
(ao, c^ii Q^j 0^3) € N 4 such that ao + a\ + a 2 + «3 ^ we have (3.20 1. In order to obtain (3.19 ), 
we just have to notice that 



E(n) = F(n) - F(n - 1) 

(n + l)(n + 2)(n + 3) 



3.2.3 A related overdetermined system 



Let us now study the equation (3.16 1, where the data is (ip 1 , ip 2 , (p 3 , ip A ) G C°°(Qo) 4 an d the 
unknown is (z 1 ,^ 2 ,^) € C°°((5o) 3 - 

Let us explain the idea behind the reasoning we are going to develop in this subsection, 



which comes from [19, Section 2.3.8, p. 148]. Equation (3.16) is "analytically" overdetermined, 
since we have more equations (4) than unknowns (3). However, if we see ( 3.16| as a linear 



system of algebraic unknowns (the unknowns being z 1 , z 2 , it and their derivatives) the system 
is now "algebraically" undcrdctcrmincd: We have 4 equations and 19 unknowns. But it is easy 



to obtain as many new equations as we want: It suffices to differentiate (3.161 enough times. 



Some new "algebraic unknowns" (the derivatives of z 1 ,z 2 ,tt) appear, but since the system was 
analytically overdetermined, one can hope that they are not too many new unknowns appearing. 
Notably, one can hope that, after differentiating a sufficient number of times, we obtain more 
equations than algebraic unknowns. We would then deduce Assumption A2 by inverting in 
some sense this well-posed linear system (this will be explained in detail later). 

We first eliminate n in our equation ( 3.16| . To reach this goal, in (3.16 1, we apply 83 to the 
first and second lines, and use the third line. We obtain the following equations: 



t2 ,1 



2/3^2 + 



(vl 



vl)4 



tz 1 
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y 2 zl 3 



-yl4 



-22,-22 
-2/3*1 +Vl z 3 



Azl 



y\4 



y\ z l 



y *i3 



l 3 z 2 - 
1 z 33 



—2 2 1 1—2 — 3\ 2 
V z 23 + (V2 ~ 2/3>3 



2 S 

^3 ~ ¥21 
2 4 

z 2 = f ■ 



"3t 



Az 2 3 



(3.21) 



The first and second equation of (3.21) contain derivatives of z 1 and z 2 up to order 3 and the 



third equation derivatives up to order 1. We would like to have the same derivatives of maximal 
order appearing in the three equations in order to be sure that the derivatives of maximal order 
appearing in the first and second equation might also appear in the third one. Hence we are 
going to differentiate the last equation 2 more times than the others. If we count the maximum 
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number of derivatives of z 1 and z 2 we create by differentiating n times the first and second 
equation and n + 2 times the third one, we obtain 



H(n) = 2F(n + 3) 



(n + 4)(n + 5)(n + 6)(n + 7) 
12 



different derivatives. The number G(n) of equations we obtain is then 

G(n) = 2F(n) + F(n + 2) 

_ (3 + n)(4 + n)(34 + 17n + 3n 2 ) 
~ 24 ' 



(3.22) 



(3.23) 



From (3.22) and ( 3.23| , one sees that G(n) — H(n) is increasing with respect to n and that 

G(18) - #(18) = -44 < and G(19) - (19) = 460 > 0. 

Hence, in order to have more equations than unknowns and as few equations as possible, we 
choose n = 19. We have G(19) = 30360 equations and H(19) = 29900 unknowns. We can see 
this system of 30360 partial differential equations as a linear system of the form 

L (t,x)Z = <Z>, 

where L € G°°(Qo; -^30360x 29900 (K)), Z € M 29900 (Z contains the derivatives of z 1 and z 2 up 
to the order 19) and $ G JJ 3 °360 ^ con tains the derivatives of ip up to the order 19). We order 
the 29900 lines of Z so that 



'Z 1 




z 2 


= 4, 


z 3 


= 4, 




' z* 


= 4, 


z 5 


= 4, 




= 4- 



Assumption _4 2 can then be written in the following form: Prove the existence of a nonempty 
open subset ujo of C 2 and of a map N £ G°°(Qo; -^6x3036o(K)) (N is the algebraic version 
of the linear partial differential operator Af introduced in Subsection |3.2| every linear partial 
differential operator can be alternatively considered as a matrix acting on the derivatives of the 
arguments of the operator) such that 



N(t, x)L (t, x)Z = (Z 1 , Z 2 , Z 3 , Z 4 , Z 5 , Z 6 ), V(t, x) € Q , VZ e 



d 29900 



(3.24) 



Since the size of the matrix Lo(t,x) is very large, it is impossible to find some N verifying 
System (3.241 by hand and we will have to do computations on a computer. Notably, it would 
be more convenient to make L be a sparse matrix in order to use relevant tools adapted to 
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the study of big sparse linear systems. This is the reason for our simple choices for a, b and c 



given in (2.16), (2.171 and (2.181 (polynomials of small order do not create to many non zero 



coefficients in Lq when they are differentiated). Using (2.1 1, (2.161, (2.171, (2.181, (2.191 and 



(2.20), System (3.211 becomes simply, in Qq, 



a{t){—Ax\ - Axix\)ez\ + a{t){-2x\x2 - 2x\)ez\ 
+a(t)(lAx\xz + 10xlx 3 )ezl + a(t){-2x\x2 - 2x\)ez\ + Aa(t)xix 2 x 3 ezl 
+a(t)(— 2x x x 3 — 2xix\x 3 )ez\ 3 + a(f)(— 2x\x 2 x 3 — 2x\x 3 )ez2 3 
+a{t)(Ax\xl + ixlxl)ezl 3 - z\ t - z\ 13 - z\ 23 - z\ 33 = (pl~ <p\, 
a(t){-2x\ - 2x\x\)ez\ + a(t)Ax!X 2 x 3 ez 3 + a{t)(-2x\ - 2xix\)ez( 
+a{t)(-Ax\x 2 - Ax\)ez\ + a(t)(lQx\x 3 + \Ax\x 3 )ezl 
J r a(t){—2x\x 3 — 2xiX2X 3 )szf 3 + a(t)(— 2x\x 2 x 3 — 2x 2 x 3 )ez 23 



+a{t){Ax\xl + Ax 2 2 x 2 3 )z 2 33 



<"3t 



-113 



-22:! 



_ 2 _ 2 _ 3 
z 333 — ^3 V21 

12 4 
-Z x - Z 2 = (fi . 



(3.25) 



Let us consider the changing of variables 



and 



ea(t) 



T-t 



(e appears when we differentiate t n- a(t) on Qq). Let M[E, S, X] be the set of polynomials 
in the variables e, s, x\, x 2 , x 3 , with real coefficients. The 30360 x 29900 entries of Lq can 
alternatively be seen as functions depending on (t,xi,x 2 ,x 3 ,e) or as elements of HL[E, S, X] 
and, from now on, we consider Lq as an element of M. 3036OX 29900(^^1 S,X]). As we will see 
after, it turns out that many of the entries of Lq are the polynomial. 

For a positive integer k, let us denote by &k the set of permutations of {1, ... , k}. To each 
cr € Sfc, we associate the matrix S a E A4k,k(^) defined by 



= i,Vie{i,...,fc}, 



Sa = o, v* € {1, . . . , k}, vj e {1, . . . , k} \ {a(i)}. 



(3.26) 



For two positive integers k and I, let us denote by Okxi the null matrix of Mkxi($t) (which is 
included in A4k X i(H[E, S, X})). The following lemma is a key step for the proof of Proposition [2] 



Lemma 1. There exist 



--(e°,S ,X°)e 
a € ©29900: 
<7 € ©30360, 
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and 

such that 



PGM 7321x7321 S, X]), 

QeM 23Q3 9x732lWE,S,X}) 

R € -M 23039x22579 (K[P, S, X]) 

a{i) =i,Vie {1,2,3,4,5,6}, 

c a _ ( P 07321,22579 \ 

SM -\q r ) ' 

the rank of P(£°) is 7321. 



(3.27) 
(3.28) 
(3.29) 



Let us assume for the moment that this lemma holds and end the proof of Proposition [2j A 
consequence of Lemma [T] is the following: 

Lemma 2. There exists a nonempty open subset ujq 0/C2, T* > and e* > 0, such that 

det P(-^—,ea(t),x) ji 0, VT e (0,T*], Vt € [T/2,T), Ve € (0,e ], Vz e w . (3.30) 



Proof of Lemma [2} Let us first point out that det P e R[E,S,X] and, by ( |3.29| , this 
polynomial is not the polynomial. Hence there exist a nonnegative integer m and a polynomial 
P € RLE, S, .XI such that 



det P(E, S, X) = S m P{E, S, X), 
P{E, 0, X) £ R[E, X] is not the polynomial. 



By (3.32 1, there exist S > 0, C > and a nonempty open subset w of C2 such that 

\P(e,0,x)\ > 2(5, Ve e [C",+oo), Vx G oj . 



(3.31) 
(3.32) 

(3.33) 



By the mean value theorem, there exist a positive integer I and a positive real number C* such 
that 



|P(e,s,x) -P(e,0,ac)| ^C*\s\ (\e\ l + \s\ l + lj , Ve <= E, VseR, Vi£w . 
By ( 2.16| , there exists e* such that 



|a(t)|((T-t)-'+ £ *VWl' + l) < ^ VTe (0,2/C], Vie [T/2.T) 



(3.34) 



(3.35) 



From (3.33), ( |3.34 I and (3.35), we obtain that 
VTe (0,2/C], Vie [T/2,T), Vee (0,e*], VaGwb, 

which concludes the proof of Lemma [2] 
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Let us now go back to the proof of Proposition [2j For every positive integer I, we denote by 
Id; the identity map of R l . By ( |3.30| ), there exists U G C°°(Qo; ■M 73 2i X 732i(K)) such that 



U(t,x)P(t,x,s) = M7321, G u . 
Let U € C 00 (Q ;X732ix3036o(IR)) be defined by 

U(t, x) := (U(t, x) 67321,23039) , V.t e w . 



(3.36) 



(3.37) 



From (3.28), (3.361 and (3.37), one has 



U(t,x)SaL (t,x) = (W7321 07321,22579) ^o- 1 , Va: G u) . (3.38) 
Let K G A^ 6 ,732i(K) be defined by 

K := (Ida 6! 73i 5 ) • (3.39) 
From ( |3.27[ ), ( |3.38| and ( |3.39[ ), one has 

KU(t,x)S a L (t, x) = (ld 6 6i 29894) S" 1 = (ld 6 6 ,29894) , Vx G LO Q , 

which shows that (3.241 holds with N(t,x) := KU(t,x)Sa, and ends the proof of Proposition [2] 



To finish the proof of Proposition [2] it is now enough to prove Lemma ^ 
Proof of Lemma [TJ The fact that the dependence of y and its derivatives in the time 
variable is quite complicated (it is both exponential and fractional) compared to the dependence 
in the space variable (which is polynomial) is problematic, because it is not very convenient to 
use for computations on a computer. In the previous proof we have seen det P as a polynomial 
in s = ea(t), e — 7^—^ (which corresponds to terms appearing when we differentiate t *-> a(t)) 
and x. Assume that we fix e = 0: This is equivalent to do as if" the derivatives of a were all 
identically the null function, i.e. to do as if the function t 1— > a(t) were replaced by a constant 
function, which is simpler than our original function a. We will then impose e° = for our 
computations. Let us set £° := (e°,s°,x°) with e° = 0, s° = 1 and x° = (1.1 , 1.2, 1.3). 

First of all, let us prove that a decomposition under the form of (3.281 exists at least at 
point £°. We present in the Appendix [a] how we computed the matrix 

^0 := L (£°) ^ -^30360x 29900 W 

thanks to a C ++ program. 

From now on we assume that we have matrix L[j at our disposal and we are going to explain 
how to exploit it to obtain Lemma [T] 

We begin with reordering the columns so that the null columns of are moved to the last 
columns. One verifies for example thanks to Matlab that there are exactly 140 such columns. 
There exist a g 629900 an d N° G A430360x2976o(KO such that 



18 



L° S^=(N° °30360xl4o) 



(3.40) 



One problem is that it could happen that some columns of Lq are equal to but the correspond- 
ing columns of Lq are not identically null. However, we check that it is not the case (thanks to 
the evaluation function described in Appendix |A| . 

Let us recall that our goal is to extract a well-chosen submatrix of Lq which is of maximal 
rank. A reasonable hope would have been that the matrix A^ (of size 30360 x 29760) itself is of 
maximal rank 29760 (we would then have obtained something similar to Lemma [I] by choosing 
some squared extracted matrix of maximal rank P° of N°, which is always possible, the matrix 
P° would then have been of size 29760 x 29760 and the non-selected lines would be permuted 
to obtain matrices Q° of size 600 x 29760 and R° of size 600 x 140) . However it turns out to be 
false, as we will see later on. 

Since computing the rank of Lq on a computer is too long because of its size, we introduce 
the notion of structural rank. 

Definition 4. Let A 6 AL„. m (lR) and B £ A^ n , m (R). We say that A and B are structurally 
equivalent if the following property is verified: 

Atj < > Bij = 0. 

This is an equivalence relation on Ai n m (M.), and we will call Cl(A) the equivalence class of A. 
The structural rank of a matrix A (denoted sprank(^4) in the following) is the maximal rank 
of the elements ofCl(A). Equivalently, if we fill randomly the nonzero coefficients of A with 
probability 1 the rank of A is the structural rank. 

One sees that the structural rank does not depend on the coefficients of the matrix but only 
on the distribution of the zeros in the matrix and is never less than the rank. The advantage of 
the structural rank is that it can be computed fast (in a couple of seconds in our case) , especially 
on sparse matrices. It corresponds to the function sprank in Matlab. 

Computing the structural rank of N° thanks to Matlab we find that 

sprank(A°) = 28654 < 29760, 

hence there is no hope that the rank of N° is maximal. 

To extract a submatrix of A which is of maximal rank, we can, for example, begin with 
extracting a submatrix of P° which is of maximal structural rank, and verify that it is of 
maximal rank too. The right way to do this is to explore more carefully how the structural rank 
is computed. In fact the key point is the existence of a decomposition in block triangular form 
(which is related to the Dulmage-Mendelsohn decomposition for the bipartite graph associated 
to any matrix, see [14] and [26 ) of a matrix. 
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Proposition 3. Let A be a matrix. Then one can permute the columns and the lines of A to 
obtain a matrix of the following form: 



(An A 12 


A 13 


A 14 \ 





A 23 


A 2i 








A 3i 


V 





A 44 J 



(3.41) 



where: 



1. (An, A12) is the underdetermined part of the matrix, it always has more columns than 
rows. 

2. (j4i3,j4i4) is the overdetermined part of the matrix, it always has more rows than columns. 

3. A12, A23, A.34 are square matrices with nonzero diagonals (notably these matrices are of 
maximal structural rank) 

4- A23 is the well-determined part of the matrix (if the matrix is square and non-singular, it 
is the entire matrix). 

Moreover, one can permute rows and columns so that A 23 is also block triangular. The decom- 
position obtained is called the block triangular form of matrix A. The structural rank of A is 
given by the sum of the structural ranks 0/ A12, A23, ^34- 



The block triangular form (3.41) (also called the coarse decomposition) of the matrix is in 
fact given by the dmperm function in Matlab, which also gives the permutation that makes the 
matrix be in the form of (3.41 ). It also gives the block form of the well-determined part (which 
is called the fine decomposition) . One easily understands how to obtain a matrix of the form of 
(3.28) thanks to this decomposition: One can permute the blocks to obtain 



(A 

A. 
A 



34 








\ 


44 











14 


All 


A 12 


A 13 


24 








A 23 J 



(3.42) 



from which we easily deduce decomposition (3.28). 

To simplify the computations, we are not going to apply this block triangular decomposition 
directly to Ljj but to a well-chosen submatrix Lq. First of all, we go back to Lo and select some 
equations and unknowns: There exist er° £ ©29900, ^ € 63036O, Q° S -^16623 x 14630 (M) and 
R° £ M 16623 x 15270 (M) such that (see ( |3.42| ) 



0l3737x 15270 

R° , 



(3.43) 



20 



where LP corresponds to the equations we obtain by differentiating the two first equations 
of ( |3.21[ ) 15 times and the last equation 17 times, so that LP is of size (C?(15), H (15)) = 
(13737, 14630) (here there are more unknowns than equations but we will see that this will not 
be a problem). 



Let us call 
and let us call 



L° :=£°«°). 

Q° ■■= Q(£°), 

R° := R(e), 
S° := 



One has 







13737x15270 



Thanks to Matlab, we find the Dulmage-Mendelsohn decomposition of L[j and observe that 



there exists some permutations matrices a 1 and a 1 such that (see (3.42 1 ) 




09050x5578 \ 

BP )' 

with Q° G A^4687x905o(1K)j R° S -^4687x5578(1^)1 an d L is of maximal structural rank and 
square of size 9050 x 9050 (it corresponds to the block A34 in the triangular block decomposition) . 

Applying the Dulmage-Mendelsohn algorithm now on L 0l we can write L in an block upper 
triangular form with 352 diagonal blocks (this is the fine decomposition), the first 351 of them 
being of "small" size and the latter one being of size 7321. Let us call L ^ ^ (with € 
[|1;352|] 2 ) the blocks of T? . 

Using this decomposition, one can see (using Matlab) that L is not of maximal rank. 
However, by computing the rank of the block L Q352 352 thanks to Matlab, one sees that it is of 
maximal rank 7321. Moreover, we verify that 

the columns corresponding to the unknowns z-y , z\ , z% , z-y , z 2 , z\ 

appear in this block, (3.44) 

by looking carefully on Matlab where the columns corresponding to these unknowns have been 
moved under the action of the permutation matrices S a o and S a i . More precisely, z\ corresponds 

to the i-th column of Lq^ 352 352 ^ and zf to the (3632 + i)-th column of LP 352,352- 
Let us call 

p o ._ 7 ° 

r -~ -^0352,352' 
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There exist a G ©29900, cr € 630360, Q° € -M 2 3039x732i(K[£, S, X]) and i?° G -M 2 3039x22579(K[£, S, X] 
such that 

&^=(£ ° 7321 ; 22579 ) < (3 - 45) 

the rank of the first block P° being maximal. The distribution of the nonzero coefficients of P° 
is given in Figure 2. 



nnz = 151219 



1000 - 



2000 



3000 



5- 4000 



5000 



6000 



7000 - 




1000 2000 3000 4000 5000 6000 7000 
unknowns 



Figure 2: Distribution of the nonzero coefficients of P°. 



Now, we come back to the matrix Lq and we consider the following matrices: 

M° := S S L°S„, 



where S& and S a are introduced in (3.45 1 



Let us call the set of all coefficients of Lq (considered as a polynomial in the variables 
E, S, X) that are not identically zero, and 0° the set of all coefficients of Lq that are not equal 
to 0. Clearly 0° C & (in fact thanks to Matlab one can see that 0° is much smaller than O), 
moreover O \ O corresponds to the nonzero coefficients of the matrix that become identically 
null when we change a(t) into the function identically equal to f and apply it at point £ . What 
could happen is that M° is not of block triangular form as in (3.45) (the null block of the 
matrix M° may contain some elements of \ 0°). However, since the only important thing is 
the location of the elements of \ O and not their value, one can verify easily on Matlab that 
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the coefficients of \ 0° do not influence the block form (3.451 (by looking where the elemnts 
of \ 0° are moved under the action of and S a ), i.e. the permutations S& and S a also give 
a decomposition under the form of (3.451 for L a : There exists 

P eX7321x732l(M[£;,5,X]), 
QeM 23 039x7321 

and 

R G A^23039x22579(R[-£', S, X]) 

such that 




(3.46) 



with the rela tions P(g°) = P°, Q(£°) = Q°, P(£°) = P°. 

Property ( |3.28 1 follows then directly from (3.431 and ( |3.46[ ), (3.291 is a direct consequence of 
the above construction, and (3.271 can be easily deduced by permuting some lines and columns 
of So and S a (thanks to Property (3.44 1). Finally, Lemma [I] holds. ■ 

Consequently Proposition [2] holds. Moreover, one observes that the linear partial differential 
operator A^o that we have just created so that (3.11 1 holds is exactly V* , (where V is the differen- 
tial version of the matrix P seen as a partial differential operator acting on (z\, z\, z\, z\, 2§> )) 



and is of order 17 (because we have differentiated the equations of (3.161 16 times, and (3.21 1 
was obtained by differentiating System ( |3.16| 1 time). Hence, the corresponding operator Ai 



in equality (3.5) is also of order 17. 



This concludes the proof of Proposition [2] 



3.3 Controllability of the linear control system (3.6) 



In this subsection we prove some technical lemmas that imply the null controllability of ( |3.6| 
with controls which are derivatives of smooth enough functions having a small support. This is 
needed to ensure that the controls are in the image of B (this is exactly Assumption .Ai) and 
to take into account Remark [3] 

The first lemma is the following one, (remind that y is one of the trajectories constructed in 
Section [2]). It consists in a Carleman estimate with curl observation. 

Lemma 3. Let 9 : f2 — > [0, +oo) be a lower semi- continuous function which is not identically 
and let r G (0, 1). There exists C\ > such that, for every K\ > C\, every v K\(l/r— 1), there 
exists s° such that for every e £ (0, e°), there exists C > such that, for every g G L 2 {(Q, T) xfi) 3 
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and for every solution z £ L 2 ((0, T), H) fl L°°((0, T), V) of the adjoint of the linearized Navier- 
Stokes system 

-z t - Az- (y ■ V"> - (z ■ V)y + Vtt = g inQ, 

V • z = in Q, 



(3.47) 



z = on [0, T] x dn, 



one has 

\\e^^z\\l 2{{T/2T)Hl{n}3) + \\z(T/2, -)\\ 2 mnr 



(3.48) 



Be ( T -') 5 |V A z\ 2 + 

(T/2,T)xQ J(T/2,T)xQ 

Proof of Lemma [3j In this proof, our system, which is initially defined on (0,T), will only 
be considered on the interval of time (T/2,T). In fact, in what follows (see, in particular, the 
proof of Proposition |4]) , we will not act on the system on the interval (0,T/2), hence we only 
need a Carleman estimate on (T/2, T). For our proof, we need to use the particular form of our 
y in time, notably that (see ( |2.16[ )), 



\y(t,x)\ + \Vy(t,x)\ < Cee^-' > 5 , V(t,x) e (T/2.T) x Q. 



(3.49) 



Without loss of generality, we may assume that there exists nonempty open subset w* of f2 such 
that 6 — l u *. Let us now give some extra notations. Let rf E C 2 (£l) such that |V?7 | > in 
D, \ui* and 77 — on dil. For the existence of 77°, see |17l Lemma 1.1, p. 4]. Let us call 

e m||Vllo 
a(t, x) := - 



gAClOH^IU+W)) e A(10|| J ,°|U+r ; (x)) 



{t-T/2) 5 {T -t) 5 



{t-T/2f{T-tf 



and 



We call 



a* (t) :— maxa(f, x). 



Q /2 := (T/2,T)xfi. 



Using (3.491 and |20l Proposition 3.1, p. 6] on the adjoint system (3.47) (where we see the 



first and zero order terms of this equation as a second member, because Proposition 3.1 of 
concerns only the Stokes system), one has, for A > C and s > C, 



S 3 A 4 



e -2«a-2.a*£3| V A z |2 + ^2 / e _2 a a-2«a*^| y(y A ^|2 
— 2sa — 2sa* /-3 



^ c I s 3 A 4 / 

'(T/2,T)xw* 



g -2 S a* e (T - t ) 6 (M : 



|Vz| 2 )+ / 



e- 2 -*| 9 | 2 
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In fact, looking carefully at the proof of Proposition 3.1 of [50], one remarks (just by changing 
the weight p(t) := e~ sa by p{t) := e~ asa where a is a parameter that can be chosen as large 
as we wish), that the previous inequality can be improved in the following way, as soon as s is 
large enough, for every a > (the constant C depending on a): 

s 3 A 4 f e- 2sa - 2asa * A z\ 2 + s\ 2 [ e- 2sa - 2asa * Az)\ 2 

<C(s 3 A 4 f e- 2sa - 2asa * f\V A z\ 2 + s f e- 2asa "e~^^(\z\ 2 

V J(T/2,T)xlu* JQ /2 



+ |Vz| 2 )+ / e- 2asa '\g\ 2 ) . (3.50) 
Jo.,-, I 



As usual, we now change our weights so that they do not vanish at time t — T/2. Let us call 
I : [T/2,T] -» R defined by l(t) = T 2 /16 on [T/2,3T/4] and l(t) = (t-T/2)(T-t) on [3T/4,T]. 
Let us define 

e 12A||r,°|U _ e A(10||Vl|oo+»7 (*)) g A(10| | V ° I loo +V° (*)) 

M,x):=- -,7(*,*):= i8(t) , 

(3*(t) :— max/3(i, x), 7*(i) := max7(i, x). 

Clearly, the functions a and /3 coincide on [3T/4, T], as well as the functions £ and 7. Using 
classical arguments (see for example the proof of Lemma 1 of [IS]), one easily obtains (one can 
forget the weight s and A in front of our terms from now on) the existence of C (depending now 
on s, A, which are assumed to be large enough, and a) such that 

e -2(l+a) S /3*| VAz |2 + f e -2(l+a) 5 /J*|V(VA*)| 2 



<C / e- 2a ^V 3 |VAz| 2 +e / e- 2as/3 V^^(|z| 2 + |V2| 2 ) 

\J(T/2,T)xu* Jq /2 



e- 2as/r |ff| 2 ■ (3.51) 

Q/2 J 

One remarks that, since V • z — in Q and z — on (0, T) x Sf2, one has 

C || e -2( W v A z|| 2 2(Q/2)3 ^ ||e - 2( W Vz||| 2(Q/2)0 . (3.52) 
Using Poincare's inequality, we also have 

C || e -2( 1+a)srv ^2 )g ^ ||e -2( Wz ||2 2( ^ )3 _ (353) 
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Putting this into (3.511, one obtains 



'Q/2 



e -2(l+a).0'| Vj8 |2 + / e -2(l+a).r| z | 



Q/2 



^ C 



-2as/3' 



(T/2,T)xw* 



7* 3 |V Az\ 2 +£ 



~2asB' 



e <^(|z| 2 + |Vz| 2 ) 



/2 



e- 2asfs "\g\ 2 



(3.54) 



Let us define, for a > 1, 



K := 2 6 (1 + a + y/a)s- 



^(io||»? ||oo+min?7 (a;)) 



x£fl 



K x := 2 6 (a - y[a)s 



A(io|| t , |U+min?7 (a;)) 

y5 



One has the existence of C* > (depending on a > 1, A >> 1 and s >> 1) such that 

e - 2os/3 * (t) 7* 3 (*) < C*e(^5)ff, Vi € (T/2.T). 
Moreover, there exists C 1 > (depending on a > 1, A >> 1 and s >> 1) such that 

Kg ^ 2(l + a) 3 /3* 

e c? 3 ^ Ce < 5 <'> , Vt€(T/2,T). 
Fixing s, A and making a — > +oo, one easily sees that Kq/K\ = 1+a+ ^ r 



(3.55) 
(3.56) 

(3.57) 

(3.58) 



X° — > 1 + so that for every 

a — \/a 17 

r € (0,1), we have for a large enough, < K\/r. For e > small enough and for v large 
enough {y ^ Kq — Ki), one can absorb the undesired terms £ J. 



-2asB* 



Q/2 



e (\z\ 2 + \V z\ 



from the right-hand side of (3.541. Then using some classical energy estimates together with 
( |3.57[ ) and fl3.58| ), one obtains ( |3.48"] ). ■ 

From now on, we set 

p r (t) := e-( T -*> 5 , /)i(t) := e^^. 

Let us now derive from this Carleman inequality a result of null-controllability with controls 
which are derivatives of smooth functions. Let l Wo : M 3 — s- [0, 1] be a function of class C°° which 
is not identically equal to and having a support included in wo, where ojq was introduced in 
Lemma [2] We apply Lemma [3] with 9 = 1 WQ . One has the following proposition. 
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Proposition 4. With the notations of Lemma^ let f G L 2 {Q) 3 be such that p r 1 / 2 f € L 2 {Q) 3 
and let us consider the following linearized Navier-Stokes control system 

y t - Ay + (y ■ V)y + (y ■ V)y + Vp = / + V A ((V A u)l Uo ) in Q, 

V ■ y = inQ, (3.59) 

y = on [0,T] x ,90, 



where the control is u. Then, for every y 6 V, there exists a solution (y,p,u) of (3.59), such 
that y(0, •) = y° and for every K\ verifying < K\ < K\, 



e 2 ( T- t) 5 (va«)1 Uo G L 2 ({T/2,T) 1 H 53 (n) 3 )r\H 27 ({T/2 1 T),H- 1 (nf), 
e^Wy e L 2 ((T/2,T),H 2 {n) 3 )nL 00 ((T/2,T),H 1 {n) 3 ). 



(3.60) 
(3.61) 

Proof of Proposition^ For what follows, we only control on the interval of time (T/2, T), 
i.e. we set it = on (0, T/2) and let the corresponding solution (y,p) of (3.59 ) on (0, T/2) evolve 
naturally until time T/2. Let y T / 2 = y(T/2 : ■). We are now going to define (y,p, u) on (T/2, T). 

Let us call X m := L 2 ((T/2, T), H 2rn (fl) 3 ) n H m ((T/2, T), L 2 (n) 3 ), which is an Hilbert space 
for the usual scalar product, and X' m its dual with pivot space L 2 (Q/ 2 ) 3 - 

Let £ : C°°(Q/2) 3 — > C°°(Q/2) 3 be the linearized time-dependant Navier-Stokes operator 

Cy = y t -Ay + {y- W)y + (y ■ V)y 

and C* : C oa (Q/ 2 ) 3 -> C°°(Q/2) 3 be its adjoint. Let P be the set of (w,h) such that 

w E L 2 ((T/2,T),H 2 (n) 3 )r\H 1 ((T/2,T),L 2 (n) 3 )r\L 0O ((T/2,T),H 1 (n) 3 ), 

h&L 2 {{T/2,T),H l {Q)) and / h{t,x)dx = 0,Vt G (T/2,T), 

V • to = in Q/2 
to = on (T/2,T) x 90, 



Let a be the following bilinear form defined on Po (according to the right-hand side of (3.48)): 
a((z, q), (w, h)) :=< y/pi(C* 'z + Vq), y/fh(C*w + Vh) >x 26 



+ I 1 U0 Pi(VAz).(VA™). 

Q/2 



From the unique continuation property for Stokes-like systems induced by (3.48), we deduce 
that a is a scalar product on P . Let P be the completion of P for the norm associated to a. 
Then, (P,a) is a Hilbert space and a is continuous and coercive on P. Let us now consider the 
linear form I defined on Pq by 



l(w, h) 



/2 



fw+ I y T ^ 2 w(T/2,.) 
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Using the Cauchy-Schwarz inequality together with (3.48), one has, for every (w,h) G Pq 



\l(w,h)\ < \\p r 1/2 /IIl2(q /2 )3||v / ^ : w|| L 2 ( q /2) 3 + ||y T/2 || L 2 (0) 3||-!«(T/2,.)|| i 2 (n) 3 

< c(\\Pr- 1/2 f\\m Q/2 r + \\y T/2 \\mnr)V<(w,h)Aw,h))- (3.62) 



Hence I can be uniquely extended as a continuous linear map from P into K. Then ( |3.62[ ) holds 
for every (w, h) G P. Applying the Lax-Milgram theorem, there exists a unique (z, q) G P 
verifying 

a{(z,q),(w,h)) = l(w,h),V(w,h) G P. (3.63) 

We then set 
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y := Pl (£*z + Vq + £(-l)*(A fc (£*z + Vg)) 



and 
where 



k=l 

S:=piVA(T W0 VAz), 



l 2 l 2 ...l k t k • 



l<ii<i 2 <...<is i <3 



Equation (3.631 implies that we have for every (w,h) G P. 

< ^(£*z + Vq),^(£*w + Vh) >x 26 + f ViVA?-VAw 



fw+ y T/2 w(T/2,.)dx. (3.65) 
Q/2 "/n 



From (3.65), one obtains that, for every g G C^{Q/2)' : 



[ fw + 


/ uw + 


J 


'Q/2 


JQo 


J 12 



as soon as (w, ft) verifies 



„T/2. 



£*w + Vh = g in Q/2, 

V • w = in Q/2, 

w = on [T/2,T] x dfl, 
{ z(T, •) = in il, 



■ X 2 



(3.66) 
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that is to say there exists a pressure p such that (y , p) is the solution by transposition of 

( % - Ay + (y ■ V)y + (y ■ V)y + Vp = f + u in Q /2 , 

V • y = in Q /2 , 

y = on [T/2, T] x dfl, 

y(T/2, •) = y T ' 2 on 17. 

Moreover, the fact that a((z, g), (z, q)) < +00 easily implies (by coming back to the definition 



(3.67) 



of the norm in X^q and using equality ( 3.64 1 ) that 

(T-t)^ Pl -^yeX' 26 . 

Using regularity results about linearized Navier-Stokes systems and a bootstrap argument, one 
obtains that, for every K\ verifying < K\ < K\, we have, with pi(t) := e ( T -*) b , 

p- 1/2 y € i 2 ((T/2, T), i/ 2 (fi) 3 ) n L°°((T/2, T), ^(fi) 3 ), (3.68) 



and (y,p), which was defined by transposition, is now the weak solution of (3.671. 

Let us now focus on u. Thanks to inequality (3.481, a((z, q), (z, q)) < +00 implies that 



p r z < 00 



(3.69) 



Q/ 



If we set u := piz, (3.691 is equivalent to 



Using regularity results for C* applied on pi 1 p^^u and a bootstrap argument, one obtains 

p^p r 1/2 u e X 27 . (3.70) 
Proposition |4] follows using (3.681 and (3.70). ■ 



3.4 Null-controllability of (3.1) 



To finish, one can gather the results of Subsection |3.2| and Subsection |3.3| in order to apply 



Proposition [T] and obtain a controllability result on (3.1 1. However, we cannot work in the C° 
setting of Proposition [T] so we need to take into account Remark [3] and to be careful concerning 
the spaces we are working with. 

Proposition 5. For every T > small enough, for every a S (0, 1), there exists ro € (0, 1) such 
that for every r £ (ro, 1), there exists C\ > such that for every K\ > C\, for every f G L 2 (Q) 

be such that e 2 *-< T -*> 5 / € L 2 (Q) 3 and for every y° G V, if 

1-r 



(3.71) 
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there exists a solution (y,p,v) of the following linearized Navier-Stokes control system 
(y t -Ay+(y- V)y + {y-V)y + Vp = f + (0, 0, l Uo v) 



V • y = 
y(0,-)=y° inn, 
y(T, •) = in n, 
y = 



in Q, 
in Q, 



on [0, T) x <90, 



(3.72) 



such that 



e m^W y e L 2 ((0,T),H 2 (il) 3 ) n L°°((0, T), iJ 1 ^) 3 ), 

»Kl 

e 2(T- t) = w e L 2 (Q). 



(3.73) 
(3.74) 



Proof of Proposition |5} We want to apply Proposition [T] First of all, we deal with 
Assumption A%. We apply Proposition [4j There exists a solution (y* ,p*,u*) of (3.591 such that 
y*(0, ■) — y° and, for every K\ verifying < K\ < K\, 



K 1 (2-l/r) 



e *(*■-'»> l UD (VAu*) e L 2 ((T/2,T),H 53 (n) 3 )nH 27 ((T/2,T),H~ 1 (n) 3 ) 



(3.75) 



e^^y* £ L 2 ((T/2,T),H 2 (n) 3 )nL°°((T/2 1 T),H 1 (n) 3 ). (3.76) 

Using well-known interpolation results (see for example [23, Section 13.2, p. 96]) and setting 
n := 27, we obtain that 

e^^T^V A u*) € H (2 "- 1)/3 (Q ) C H 17 (Q Q ). 
Let us call u> := 1oj (V A w*), which is supported in Qq. One observes that 

— Wr, 

V A w — \ wl - wf 



Hence in view of equality (3.8 1 and setting 

f l = 0, f = -w 3 , f 3 = w 2 ,f = -w 3 , f = 0, f = w\ f = w\ - wl 

one has V A w £ Im(B) and Assumption A\ holds. 

Now, we observe that Assumption A2 follows from Proposition [2] Let (y,p,v) be defined by 

{y,P,v) := -Mw, 

where M. is as in (the proof of) Proposition [2j It makes sense to apply M. to w because M is 
a partial differential operator of order 17 and w € H 17 (Q ). 
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Using the fact that the operator M. is a partial differential operator of order 17 and that 

732 li; 

the coefficients of M explodes at time t = T at rate at most e ( T -*> 5 , as it follows from the 
construction of y given in Section [2] (see in particular (2.16 1 ) and the construction of M given 
in the proof of Proposition [2] one has 



K-2 



e W^WyeL 2 (Q ), 



(3.77) 
(3.78) 



for every K 2 < Ki(2 — 1/r) — 7'521v. In order to obtain y(T, .) = 0, it is enough to have an 
exponential decay for y at time T, i.e. to impose 

Ki(2 -1/r)- 7321v > 0, 



which, with (3.711, is equivalent to 



r > 



7322 
7323' 



which can be insured since r can be arbitrarily chosen close to 1. Let a € (0, 1). We take 

7322 



r := 



7323 - a ' 



Then if r G (tq, 1), one has 



aK x <K-l- 7321^1- 



1 



By (3.801, there exists K 2 such that 



aK x <K 2 <K x ~ 7321^ 



1 — r 



K x - 732 liA 



(3.79) 



(3.80) 



(3.81) 



Finally, one can apply Proposition [T] and we set 

(y,p) ■= (y* +y,p* +P)- 



Thanks to ( |3.75| , ( |3.76[ ), ( |3.77| and ( |3.78[ ), one has 

e^c^u e £ 2 (Q ), 
e^r^y e L 2 (Q ). 



Then, using usual regularity results for the linearized Navier-Stokes operators on e 2 ( T - t)5 y 
now considered on the entire time interval (0,T) (y € L 2 ((0,T/2) x f2), hence we also have 

e 2 < r -') 5 y £ L 2 ((0,T/2) x fi)), we obtain 

eic&y g L 2 ((0, T), iJ 2 (ft) 3 ) n L°°((Q, T),H\n) 3 ), 
as soon as y° E V. The proof of Proposition [5] is completed. ■ 
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4 Proof of Theorem [T] 

To conclude, we are going to apply an inverse mapping theorem to go back to the nonlinear 
system, which is the following (see [5] Chapter 2, Section 2.3]): 

Proposition 6. Let E and F be two Banach spaces. Let eo G E and A : E — > F which is 
of class C l in a neighborhood of eo- Assume that dA(eo) G C C {E,F) is a surjective operator. 
Then there exist r\ > and C > such that for every g G F verifying \ \g — _A(eo)|| < rj, there 
exists e G E such that 

1. A(e) = g, 

2. ||e-e |U<C|| 5 -^l(eo)|| F . 

We are going to use the same techniques as in [T2]. Let a G (0, 1), and let us set consider 
a r G (ro, 1) where ro verifies (3.791. We apply Proposition [6] with E and F defined in the 
following way. We let E be the space of the functions 

(y,p,v) G L 2 {Qf x L 2 (Q) x L 2 {Q) 

such that 

1. e^^y G L DO ((T/2,T),ff 1 ((l)) 3 nI 2 ((T/2,T) ) ff 2 (n) 3 ), 

2. V • y = in Q, 

3. VpG L 2 (Q), 

4. e^-^v G i 2 (Q) 3 and v = on Q \ [T/2,T] x w , 

5 . ( yt _ Ay + (y • V)y + (y ■ V)y + V P - (0, 0, «)) G £ 2 (Q) 3 , 

6. »(o,-)ey, 

equipped with the following norm which makes it a Banach space: 

\\(v,P,v)\\e = ||e 2 ( T - t ) 5 2/||L~((o,T),ffi(n)3)nL2(( ,T),H2 (!:2)) 3 + |b|| L 2(Q) 
+ || e i#=Vi;|| i2(Q)3 +||eM^) r (y t -A+(y-V)y+(2/-V)y+Vp-(0, 1 «))|| L » ( Q ) s+||v(0 > -)||Hi(n)». 
We let F be the space of the functions (h, y°) G L 2 (Q) 3 x V such that 

e 2(T-t)5 /j g L 2 (Q) 3 , 

equipped with the following scalar product which makes it a Hilbert space: 

((ft, j, ) |(M°)) = (e^^h\e^Wk) L2(Q)3 + (y \z°) HHQ) s. 
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We define 

A(y,p, v) = (y t -Ay + (y V)y + (y ■ V)y + (y • V)y + Vp - (0, 0, v),y(0, ■)). 
To apply the previous inverse mapping theorem, we first show the following lemma. 
Lemma 4. The map A has its image included in F and is of class C 1 on E. 

Proof of Lemma |4j 

We see that A = B + C with 

B(y,p, v) := (y t - Ay + (y ■ V)y + (y ■ V)y + V P - (0, 0, »), y(0, ■)). 

and 

C(y,p,v) := ((y V)y,0). 

Thanks to the construction of E and F we have B : E —> F and B is continuous, so, since B 
is linear, B is of class C . The map C is a quadratic form, hence to prove that it maps E into 
F and is of class C , it is sufficient to prove that it is continuous, i.e. to prove that 

\\e^^(y-V)y\\ L 2 (Qor ^C\\(y,p,v)\\%. (4.1) 
We choose r and a so that 



r 



< 2aK l . (4.2) 



(One can take for example a = 3/4 and then r £ (0, 1) close enough to 1.) Let us call 



y(t,x) := e^^-^y. (4.3) 



This definition of y and inequality (4.2) imply that 

l|y|U~((o,T),ffi(n) 3 ) sC C\\(y,p,v)\\ E , 

which gives that 

\\V(y)\\L~((o,T),mn)°)<C\\(y,p,v)\\ E . (4.4) 

We also have 

l|y||i=((o,r),H 2 (n)3) < CIK^P^)!!^- (4-5) 



A classical Sobolev inequality in dimension 3 together with (4.5 1 imply that 

l|y||L2((o,T),L~(ri)3) «S C\\(y,p,v)\\ E . (4.6) 
Direct computations imply that 

\\{y ■ V)y|| L 2 ((0iT)ii 2 (o))3 < ||V(y)|| L oo ((0:T)iL 2 (n)9) ||2/|| L 2 ((0:T)iL oo (n)3) . (4.7) 



From ( |4.4[ ), ( |4.6| and ( |4.7| , we obtain 

IKS ' V)y|| L 2 ((0iT)iL 2 (a)3) < C||(?/,p,w)|||, 



which, together with (4.3|, gives (4.1). 
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We now consider the element eo = (0,0,0) and we compute 

<U(eo)(tf, q, v)=y t -Ay+(y- V)y + (y ■ V)y + V P - (0, 0, v). 

The application of Proposition[5]we gave in the previous subsection implies that this application 
is surjective. Hence, taking g — (0, y°) and applying Proposition [6j Theorem [l] easily follows 
(notably because the trajectory y can be chosen as small as we want since e can be arbitrarily 
chosen small). 



A Appendix: Creation of the matrix Lq 



In this appendix we explain how the matrix Lq at point £° which represents all the differ- 
up to the order 19 was created. 



entiated equations of System (3.25 

C++ 



The program is written in 

using the library uBLAS which is well-adapted to the manipulation of sparse matrices. 
It is a parallel openMP algorithm, using 8 cores. We are not going to give all the technical 
details but just explain rapidly the spirit of the algorithm. To simplify, we will assume that the 
following "black boxes" (that had to be created) are at our disposal: 

1. An evaluation function ep which evaluates a polynomial (represented by a vector) at 
£°. This evaluation function can be created so that it can verify that £° is not a root of 
the polynomial P(0, ., .). (one just has to see if the evaluation is equal to whereas the 
polynomial has nonzero coefficients). 

2. A derivation function deqex which differentiates an equation of level m with respect to 
Xi,x 2 ,%3 or t. 

A partial differential equation which is a derivative of order m of some of the equations of 



(3.25) will be represented in a matricial form in the following way: We know that there are at 
most F(n + 3) derivatives appearing, and we observe that the coefficients are polynomials in 
(xx, x 2 , x-i) of an order less than 4 (it is a vector space of dimension 35). Hence the equation is 
represented by a matrix with F(n + 3) lines and 35 columns, where on each line one can find 
the coefficient of the partial derivatives of z 1 (or z 2 appearing) corresponding to the number of 



this line, thanks to a natural bijection between N 4 and N. Since we have 3 equations in (3.251 
and 2 unknowns (z 1 and z 2 ), one can write the matrix M in the following way: 



Mi 


Bi 




D 2 


M3 


B 3 



For i = 1,2,3, Ai represents the derivatives of z 1 appearing in the derivatives of the i-th equation 
of (3.25 1 and Bi those of z 2 . Hence, we can compute these At and Bi separately and then gather 



them to obtain Lq. 



The algorithm is the following. We explain it for the first equation of (3.25) and for the 



unknown z (i.e. for A±, but it is the same for the other matrices). 
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1. We create a matrix e that represents the equation. We use ep to fill the line of Lq 
corresponding to the equation in a .txt file under the form i j Ai(i,j). We create a matrix 
h which is empty for the moment. In fact in e we will keep the equations of level m — 1 
and in h we will fill the equations of level m. 

2. We create a "for" loop on m which will represent the level of equations we are creating. 
The integer m goes from 1 to 19 since we differentiate 19 times at most. 

3. We create a second "for" loop in the interior of the first loop on a number n which represents 
one of the equations of level m. Thanks to the definition of the function F given in 
Subsection 3.2.2| we have F(m — 1) + 1 ^ n ^ F(m). If m = 1, then n goes from 



F(0) + 1 = 2 to F(l) = 5 (n represents d\, 82, d% or dt). If m — 2, then n goes from 
F(l) + 1 = 6 to F(2) = 15 (n represents d\ x , df 2 , df 3 , Sf 2 , 9| 3 , 3f 4 , 5| 3 , c& or 9 t 2 t ), 
etc. This loop is parallelized on our 8 cores. In this loop, we want to create the ri-th 
equation denoted E n , which is of level m. Hence we take a suitable equation of level m — 1 
denoted E r which is so that if we differentiate E r with respect to 1,2,3 or t, we obtain 



E n . For example, if we consider vn — 2 and we want to obtain the first equation of (3.251 



differentiated two times with respect to 1, then we consider the equation E r to be the 



first equation of (3.25) differentiated one time with respect to 1 and differentiated with 



respect to 1 to obtain E n . 

4. Once the loop on n is ended, we have in our matrix e all the equations of level m — 1 
and in h we have just created all the equations of level m. Now we just have to use our 
evaluation function ep on h to obtain the coefficients of the lines of A\ corresponding to 
the equations that are of level m, i.e. the equations numbered from F(m — 1) + 1 to F(m). 
We write these coefficients in our .txt file in the form i j Ai(i,j). 

5. We update now e, take e = h, we empty h and we can go to the following loop m + 1. 

At the end we have created a file containing the coefficients of a sparse matrix A\ of size 
(8855,14950). Using the same program with z 2 and the two others equation we obtain five 
other files representing five matrices that we gather to obtain the matrix L (^°) = Lq. Our 
matrix Lq, which represents all the equations, is of size (30360,29900) and has 651128 nonzero 
coefficients. Only 0.0717% of the coefficients are different from 0, with an average of 21.45 
nonzero coefficients on each row, which is logical since we are working with coefficients that are 
polynomials of small degree, so we do not create many terms on each line when we differentiate 
the equations. In Figure [A] one can observe how the nonzero coefficients of Lq are distributed. 



B acknowledgements* 

The authors would like to thank Sergio Guerrero for fruitful discussions concerning Propo- 
sition |4j 



35 



nnz = 651128 




Fi gure 3: Distribution of the nonzero coefficients of L\ 
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